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Abstract 

We explore the possibilities of creating radiatively stable entangled states of 
two three-level dipole-interacting atoms in a A configuration by means of laser 
biharmonic continuous driving or pulses. We propose three novel schemes for 
generation of entangled states which involve only the lower states of the A 
system, not vulnerable to radiative decay. Two of them employ coherent 
dynamics to achieve entanglement in the system, whereas the third one uses 
optical pumping, i.e., an essentially incoherent process. 
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The concept of quantum entanglement, one of the most intriguing properties of multipar- 
tite quantum systems Q, has been intensively exploited over the last decade in connection 
with quantum information processing. It has been shown that the use of entangled states 
opens new horizons in such practical fields as cryptography ||, computing ||, information 
transmission |4j], and precision measurement ||. However, all of these applications become 
possible only with a reliable source of entanglement. Traditionally, entangled particles have 
been generated in the down-conversion nonlinear process 0,[7|, but this method is in some 
cases disadvantageous due to the speedy nature of the produced particles (photons) and the 
intrinsic randomness of their appearance times. That is why efforts are now being made to 
find ways for controlled production of entangled states of less volatile massive particles || . 
During the last few years, various methods for creation of entangled states of atoms, ranging 
from continuous observation of radiative decay |],|n| to controlled cold collisions []TT|, have 



been proposed and some of them experimentally demonstrated [|12| . p~3 



Though the resonant dipole-dipole interaction (RDDI) was suggested for realization of 
entangling dynamics as early as 1995 f[4 |, it was only recently that several authors fI7^ . [l6| . |T7 



investigated this interaction in more detail as a method for entangling neutral atoms in 
optical traps (neutral atom realizations benefit from the fact that neutral atoms are less 
sensitive to stray EM fields — a major source of decoherence in ions |TB1]). While the authors 



of Ref. [15] offered qualitative arguments for realization of this idea in dipole traps [19 
the papers [I6|.[l7j considered quantitative models of creation of maximally entangled states 
of two-level atoms. Unfortunately, such entangled states of two-level atoms have short 
lifetimes due to radiative decay. Obviously, since radiative decay and the RDDI have the 
same physical nature, we cannot avoid the former while making use of the latter. In this 
paper we solve this conceptual problem by presenting new methods for creation of radiatively 
stable entanglement in a system of dipole-interacting three-level atoms. Though the model 
considered here is still far from representing the real situations (see for details of a 
possible experimental realization), it offers new insights into how the RDDI can be used to 
entangle real, multilevel atoms. 



I. THE MODEL 
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we consider here two identical three-level atoms 
The dipole transitions |1) <-> |3) and 



Extending the model described in 
in a A configuration (Fig. P fixed at a distance R 
1 2) <-> 1 3) of both atoms are driven by two near-resonant laser fields. Taking the two limiting 
cases, we consider only two types of geometry: when the laser fields are either perpendicular 
or parallel to the radius vector R connecting the atoms (these geometries are shown in 
Fig. ^ and identified as symmetric and antisymmetric, respectively). Within the interaction 
picture and rotating wave approximation, the evolution of the system interacting with the 



laser fields is governed by the following master equation [BT 
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where the upper indices, i and j, number the atoms, the lower ones, k3 and 3k (k = 
1, 2), refer to dipole transitions of the atoms, and cr^ denotes the Heisenberg transition 
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operators from level \k) to level \l) within the zth atom. Relaxation effects in the system are 
characterized by the single-atom decay rates, 7^3 = 7^3 = 7*3 , which correspond to the 
conventional radiative decay into free space, and the photon exchange rates, 7^ 3 2 ' > = 7^3 , 
which describe collective relaxation, a well-known companion of the RDDI. The effective 
Hamiltonian 7i e ^ includes interaction with the laser field and the RDDI coupling on both 
transitions: 

n cS = h £ L 3 nf + ?fatl + Xn*%*$! + H. c.) , (2) 

i,fc=l,2 \ / 

where stands for the population operator of the level \k) in the zth atom, ^3 are the 
detunings of the laser field frequencies from the corresponding transitions \k) <-> |3) of an 
isolated atom, Q^l is the Rabi frequency of the laser field acting on the \k) <-> |3) transition of 
the zth atom, and Xk3 is the RDDI coupling strength on the \k) <-> |3) transition. Throughout 
the rest of this paper we will consider the case of wide homogeneous laser beams, so that 
the Rabi frequencies acting on the two atoms may differ in phase but not in magnitude, 

I i4 fc3 I ~~ I si fc3 I- 

(12) (21) 

Normalizing the RDDI parameters, Xk3, T k3 , and 7^.3 , by the decay rate of an isolated 
atom, 7fc3, we introduce the dimensionless parameters 

9k3 = 7fc 3 2) /7fc3 = 7fc3 1} /7fc3, fk3 = Xk3/lk3, (3) 



which are given by the following expressions [21]: 
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where e*j (i — 1, 2) is the unit vector in the direction of the dipole moment matrix element of 
the corresponding transition \k) <-> |3) of the ith atom, e R is the unit vector in the direction 
of R, and (p k3 = k k3 R is the dimensionless distance between the atoms (k k3 = io k3 jc is the 
wave number associated with the transition \k) <-> |3) of an isolated atom). Throughout 
the following discussion we will assume, for the sake of simplicity, that the dipole moments 
are real, collinear with each other, and perpendicular to the radius vector R (other dipole 
moment orientations lead to qualitatively the same results). 

In the case of two-level atoms, the simplest description of the system dynamics is 
offered by the basis of the Dicke states, which is formed by the doubly excited state, 
l^e) = |e) 2 , the ground state, \$> g ) = \g) 1 \g) 2 , and the two singly excited maximally en- 
tangled states — the symmetric, \^/ s ) = ^{\g)\ l e )2 + l e )i IflWj anc ^ the antisymmetric one, 
|^a) = ^(\9)i l e )2 ~~ l e )i 1^)2) t^ ne corresponding energy diagram is shown in Fig. ^|(a)]. For 
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the case of three-level atoms considered here it is useful to introduce simple generalizations 
of the Dicke states. The role of the ground and doubly excited Dicke states is then played 
by the three tensor product states \kk) — \k) 1 \k) 2 , k = 1, 2, 3, while the symmetric and anti- 
symmetric Dicke states now are represented by the three symmetric and three antisymmetric 
maximally entangled states \s k i) = 75(^)1 IO2 + IO2) and l a «) = vfd^i IO2 ~~ IO2)' 
k, I = 1,2,3; k < I. The corresponding energy diagram is shown in Fig. |3](b). Note 
that in both two- and three-level models the energy levels can be grouped according to 
their type of symmetry: the unentangled states \kk) = \k) 1 \k) 2 , as well as the states 
\ski) = 7^(l^)i IO2 + IO2)' can ^ e sa ^ to belong to one type of symmetry (symmetric 
with respect to the atom interchange), and the states |a^) = -^(\k) 1 \l) 2 — \k) x |/) 2 ) to an- 
other (antisymmetric with respect to the atom interchange). The transitions between these 
levels can then be classified as symmetry-preserving and symmetry-breaking, respectively. 
It is easy to show that, due to the form of the transition matrix elements, the symmetry- 
preserving transitions are sensitive only to the sum of the Rabi frequencies, ^3 +^3*, acting 
on the atoms, and the symmetry- breaking transitions only to their difference, — fijj.g . 
In the following, we also assume that the system is initially stored in the |11) state, which 



can easily be achieved by conventional optical pumping methods [22 



II. COHERENT ENTANGLING PROCESSES 

A. Resonant Raman pulses 

In our previous paper [|l(| we have shown that the maximally entangled Dicke states \^f s ) 
or \^f a ) of two two- level atoms can be efficiently populated at small interatomic distances 
simply by applying an appropriately tailored laser pulse. Assuming that initially the entire 
population of the system is concentrated in the ground state \^ g ), this pulse should be 
tuned into resonance with a transition to only one of these maximally entangled states. 
Then, by applying a pi pulse analog, a significant part of the population of the system can 
be transferred to one of these states, thereby creating entanglement in the system. We have 
also shown that the entanglement fidelity, defined as the population of the corresponding 
maximally entangled state, can be made arbitrarily close to unity as the interatomic distance 
R goes to zero. 

In this paper we propose new ways to create stable entanglement in a system of two 
three-level atoms. To be radiatively stable, the created entangled states should involve only 
the lower levels |1) and |2) of the original A system of each atom, as only these states are 
not vulnerable to radiative decay. Therefore, our goal here will be to achieve the maximum 
possible population of one of the maximally entangled states |ai 2 ) or |s 12 ) (see Fig. 0(b)). 
The most straightforward way to do this is to extend the results of the two-level model 
to the three-level one, considered here, by using resonant Raman pulses. By the latter we 
mean a sequence of two coherent pi pulses, the first of which transfers the population to 
one of the maximally entangled states involving the initial lower level of the A system and 
some quickly decaying upper lying "transit" level, while the second one transfers the entire 
population of the "transit" level to another radiatively stable lower level of the A system, 
thus removing the radiative instability of the entanglement. 
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In the considered system of two dipole-interacting three-level atoms, the role of the in- 
termediate "transit" state can be played by the above-mentioned levels |a 13 ) or |s 13 ) (one 
should not forget that the system is initially in the state 1 1 1) ) . During the first step, the 
pulse resonant, for example, with the |11) «-> |si 3 ) transition transfers the population to the 
maximally entangled state | S13) ; the second step creates the radiatively stable maximally 
entangled state |si2) by application of the symmetry-preserving pi pulse resonant with the 
I S13) <-> I S12) transition. In fact, it is the symmetry preservation rules that prevent popu- 
lation from going into the \ci12) state in a transition that is also resonant with the second 
pulse. 

For both pulses to be resonant, the parameters of the laser field should be chosen in the 
following way 

«fc3 = 0, 43 = X13/2, Ifijial < Ixisl, ( 5 ) 

where 0^3 is the phase difference between the Rabi frequencies acting on the two atoms, 
f2fc3 = ex P('i> a k3) (considering laser beams formed by traveling waves, 0^3 varies from 
zero for the symmetric geometry to ^3 for the antisymmetric one, and takes all the interme- 
diate values for other types of laser field geometry). The parameters given by (|3p correspond 
therefore to the case when both lasers are used in the symmetric geometry, the "transit" 
state is |si3), and the final radiatively stable maximally entangled state is |si2)- Other types 
of geometries and laser parameter sets can obviously be chosen when using the other inter- 
mediate state (013) and/or creating the other radiatively stable maximally entangled state 
|a 12 ) • For example, to create the |a i2 ) state, one can use the following set of parameters: 

"13 = (Pxs, "23 = 0, 8 k3 = -X13/2, < |xia|- (6) 

The phase differences a k 3 in this case correspond to one of the lasers beams being used 
in the antisymmetric geometry, and the other in the symmetric one. Note that, when using 
antisymmetric geometry at small interatomic distances (^3 <C 1), most of the laser power 
is "wasted" since only a fraction of it contributes to the corresponding transition matrix 



element 



11| Heft/h 1 0,13) = |f2^ — f^i3 |/2 = |^i3 I sin(<y9i 3 /2) <C |^i3 I, and actually induces 



transitions |11) <-> |a, 13 ) . 

While a simple estimate of the resulting fidelity of creation of the maximally entangled 
state is offered by a product of the fidelities of each step of the resonant Raman process 
(which were calculated in ]16| within the two-level atoms model), rigorous results can be 



obtained only by explicit solution of the corresponding master equation. Due to the high 
dimensionality of the master equation (|1]), this calculation is rather demanding computa- 
tionally, and was not included in the present treatment. 



B. Stimulated Raman adiabatic passage 

Another coherent method for creation of maximally entangled states is based on the 
stimulated Raman adiabatic passage (STIRAP) technique [p3|, a well-known alternative to 
Raman pulses. The STIRAP method uses adiabatic following of the system state after the 
slowly changing parameters of the laser field, which are chosen to form a so-called counter- 
intuitive pulse sequence. The STIRAP technique benefits from extremely low probabilities 
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of losing coherence due to radiative decay of the intermediate states, and has already been 
proposed for use in entanglement-related problems W% - in our case, efficient transfer of the 
population from the state |11) to the state \a 12 ) or |s 12 ) may be deterred by existence of 



several intermediate states |25) . However, as we show below, efficient transfer is still possible 
for appropriately chosen laser pulse parameters. 

To realize STIRAP in our system we need to choose the frequencies and geometries of 
the two constituent laser pulses in a way that would leave active (i.e., resonant and having 
strong transition amplitudes due to the use of the corresponding geometries, see section H) 
only two transitions in the whole system. An appropriate choice is given by 

"13 = 0, a 23 = 7T, S k3 = X13/2, max = fi < \xw\, ( 7 ) 

where f2 stands for the amplitude of the corresponding constituents of the counterintuitive 
laser pulse sequence. The condition a 23 = ir, which is very important as it prevents leakage of 
population into other levels, can be easily realized by using two laser beams in antisymmetric 
geometry, which form a standing wave with one of the nodes situated exactly in the middle of 
the vector R connecting the two atoms [I7|. In this case only two transitions, |11) <-> |s 13 ) 



and |si 3 ) <-> |ai 3 ), are active, and the adiabatic passage results in transfer of the total 
population to the radiatively stable state |ai 3 ). 

We have numerically calculated the final population (fidelity) of the state |ai 2 ) after the 
STIRAP procedure by explicit solution of the corresponding Schrodinger equation with the 
Hamiltonian given by (0). The two laser field pulses had the same Gaussian form and were 
delayed with respect to each other by their length ||23|| , and the rest of the parameters were 
given by (|7|). For determinacy, the length of the pulses was chosen to be equal to one- 
tenth of the lifetime of the excited level |3) of the original A system, t p = 0.1/(7i 3 + 7 23 ). 
The final population of the level |ai 2 ) is shown in Fig. M function of the pulses 

Rabi frequency amplitude, Q , for different values of the RDDI splitting parameter f 13 . 
As one can see from the figure, for sufficiently high RDDI splittings (i.e., for sufficiently 
small atomic separations) the fidelity first grows with increasing Q , reaching saturation at 
Q T P ~ 5, which corresponds to the adiabaticity condition on the pulse area p3| . Then, 



after some point, the final inequality in ([?[) is no longer fulfilled and the efficiency of the 
process degrades due to nonresonant excitation of other levels caused by power broadening. 
For the same reasons, the fidelity does not reach unity at any values of Qq for low RDDI 
splittings (large interatomic separations). In Fig. ^(b) we show the overall fidelity of the 
STIRAP method for optimized values of the Rabi frequency amplitude as a function of the 
interatomic distance y?i 3 . 

As we ignore relaxation processes in this model (a common practice for STIRAP sim- 
ulations), one should beware of relaxation- induced errors. However, these errors assume 
significant values only for the case of long pulses, r p > 1/(713 + 723), and low overall STI- 
RAP process fidelity, i.e., situations that are not of great concern to us here. 



III. AN INCOHERENT ENTANGLING PROCESS: OPTICAL PUMPING 

An interesting alternative to the coherent methods can be offered by optical pumping 
schemes where the stationary state of the system corresponds to one of the maximally 
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entangled states. In this situation, the population of the system is pumped into the entangled 
state after asymptotically large time periods. 

Consider the following choice of the laser field parameters: 

«fc3 = 0, 5 k3 = Xk3/2, |fife3| < \Xk3\- (8) 

Neglecting nonresonant excitation at small interatomic distances, only a few transitions 
remain resonant and have the corresponding geometry. These active transitions are shown 
in Fig. ||(a) (the upper state 1 33) is omitted in the figure, as it is only negligibly excited at 
small interatomic distances [16|]). As seen from the figure, the maximally entangled state 
1 012) is not included in the chain formed by the laser- induced transitions; however, it is still 
populated as a result of the decay of the upper-lying levels, as shown by the dotted lines in 
the same figure. As the state |ai 2 ) is stable with respect to both the laser-induced transitions 
and radiative decay, all of the population will eventually be pumped into this state. One 
should note, though, that as the interatomic distance goes to zero, the symmetry-breaking 
decay rates decrease, which leads to a corresponding increase of the required pumping time. 
If we choose another configuration that uses antisymmetric standing-wave geometry of the 
laser beams [Fig. [|(b)], 

a k3 = 7T, 5 k3 = ~Xfe3/2, |0fc3| < \Xk3\, (9) 

the increase of the pumping time is still brought on by the decrease of efficiency of the 
symmetry-breaking laser-induced transitions at small distances since the corresponding 
transfer matrix elements are proportional to Ifljjjg — }\ ~ sin (^3/2). 

Strictly speaking, the above arguments hold only in the case when the RDDI coupling 
constants \kz on different transitions are equal (possibly, up to an error on the order of 
7fc3). This condition is satisfied, for example, when the two lower levels of the original A 
system are sublevels of the same atomic level. However, even when the RDDI coupling on 
two transitions differs considerably, the present treatment is still applicable provided that 
one uses four lasers instead of two to satisfy all of the resonance conditions for transitions 
shown in Fig. |5]. In contrast to the methods presented in the previous sections, it is also very 
important to avoid high a degree of mutual coherence of the components of the biharmonic 
laser pumping, as otherwise the population of each atom will be trapped in a corresponding 



dark state 26 



To prove the foregoing arguments, we have numerically calculated the stationary states 
of the master equation ([1]) with the laser pumping parameters given by (|D and @. In order 
to disrupt trapping of the populations of the two atoms in the single-atom dark states, we 
introduced additional elastic dephasing of the lower level transition |1) <-> |2) in both atoms, 
which can be easily realized by the relative jitter of the two pumping laser frequencies. 
Assuming that this elastic dephasing is characterized by the rate r 12 , the corresponding 
relaxation superoperator, which should be plugged into the master equation ([!]), has the 
form 

fitter P = T 12 J2 {2a^pajp-paPajp-a^ajpp) 1 (10) 

where a® = — n 22 is the lower level population difference operator in the zth atom. 
For simulations we used a realistic value T 12 = O.OI7, where we again assume for simplicity 
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7 = 7i3 = 723- The results of the numerical calculations of the steady state population of 
the level |ai 2 ) for different values of laser pumping Rabi frequencies Q (in our calculations 
they are equal for all transitions and atoms, — Q, i,k — 1,2) are presented in Fig. |6| 

as a function of the interatomic distance <pi3 for the two geometries discussed. As seen 
from the graphs, the fidelity of the entanglement produced first monotonically decreases 
with increasing Rabi frequency, and then strongly degrades when the magnitude of the 
Rabi frequency approaches that of the RDDI splitting due to power-broadening-induced 
nonresonant excitation. The graphs for Q = O.OOI7 in Fig. |6|, therefore, decently represent 
the overall fidelity of the optical pumping method. For low Rabi frequency amplitudes, the 
antisymmetric geometry clearly shows better results and achieves fidelity of 0.8 at y9 13 « 1, 
which is much better than the fidelities achieved by other methods at such distances. 



IV. CONCLUSIONS 

We have considered three methods for creation of radiatively stable entanglement in 
a system of two dipole-interacting three-level atoms in a A configuration. It was shown 
that the radiatively stable maximally entangled states, \au) and | S12) , which involve only 
the lower levels of the original A systems, can be efficiently populated at small interatomic 
distances by employing coherent or incoherent methods. 

The first of the coherent methods, which employs resonant Raman pulses for trans- 
fer of population (first to a radiatively unstable maximally entangled state and then to a 
stable maximally entangled state), makes use of specific resonance conditions and symmetry- 
preservation rules. The second coherent method, which utilizes a STIRAP process, realizes 
adiabatic transfer of the population of the system into the final state coinciding with one of 
the radiatively stable maximally entangled states. The STIRAP method, however, requires 
the use of standing waves to avoid leakage of population into unentangled states. 

As a rather surprising result, we have also shown that entanglement can be deterministi- 
cally created as a result of an incoherent process [f27| , optical pumping in our case. Creating 



a laser field configuration where one of the maximally entangled states, \au) or |si 2 ), is 
not included in the chain of laser-induced transitions, we achieve high populations of that 
state at asymptotically large times due to radiative decay into that state. An important 
restriction for realization of the optical pumping method is that one has to avoid high mu- 
tual coherence of the pumping laser beams, but this restriction becomes an advantage when 
realizing the proposed schemes experimentally as it is usually easier to provide an incoherent 
pumping than a coherent one. 

For two of the proposed methods (STIRAP and optical pumping), the fidelity of the 
created approximations of the maximally entangled states was calculated, and it shows 
qualitatively the same dependence on the interatomic distance R as in the previously con- 
sidered two-level atom model fl6| . The fidelity of 0.8 (a good benchmark for Bell inequality 
violations) is achieved in all of the considered methods at interatomic separations between 
one-fifteenth and one-sixth of the wavelengths of the working transitions. 

In conclusion, we have shown that radiatively stable maximally entangled states can be 
created in a system of two dipole-interacting atoms under conditions that can be exper- 
imentally implemented, for example, in optical lattices EREDl. The general form of the 
RDDI operator also suggests that simple analogs of the proposed methods can be employed 



S 



in other physical systems, such as quantum dots in semiconductors p7| , p8| or cavity QED 
systems [^,^] (or, indeed, a combination of the latter two [30]). 
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FIGURES 




FIG. 1. The level structure of an isolated three-level atom in a A configuration. The dipole 
transitions, |1) <-> \2) and \2) <-> |3), are driven by two laser fields, which are detuned by £13 and 
(523, respectively. Dotted lines show radiative decay channels and their corresponding rates. 
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FIG. 2. Geometry of the model with directions of laser beams for the "symmetric" and "anti- 
symmetric" laser beams. 
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FIG. 3. (a) Energy levels of two dipole-interacting two-level atoms (Dicke states). The two 
maximally entangled Dicke states, |Vl/ ) and \^ s ), are split by the RDDI coupling strength \- Also 
shown in the figure are the radiative decay channels with the corresponding decay rates, (b) The 
same for two dipole-interacting A systems. Shown are the doublets, formed by symmetric and 
antisymmetric Dicke-like states, and the laser-induced transitions induced by the two components 
of the biharmonic driving. 
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FIG. 4. (a) Population of the maximally entangled state \cL12) after adiabatic passage versus 
the pulse area for different values of the RDDI parameter / = /13 = X13/713. (b) The same for the 
optimal value of the laser pulses area, £Iqt p , versus the interatomic distance 9313. In both graphs 
we assume equal decay rates of the two channels of the original A system, 713 = 723. 




FIG. 5. Energy level diagram and transitions for two variants, symmetric (a) and antisym- 
metric (b), of the optical pumping scheme. Solid lines indicate laser-induced transitions with the 
corresponding phase differences a between two atoms. Dotted lines show the significant decay 
channels with the corresponding decay rates. The negligibly populated upper state 1 33) is omitted 
in both figures. 
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FIG. 6. Population of the maximally entangled state \a\2) in the stationary solution of the 
master equation as a function of the interatomic distance 9913 , for different values of the Rabi fre- 
quencies 1 = fi, i, fe = 1, 2 and the two discussed geometries, symmetric (a) and antisymmetric 
(b). 
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